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Abstract. Carbon and nitrogen impurity transport coefficients are determined from 
gas puff experiments carried out during repeat L-mode discharges on the Mega-Amp 
Spherical Tokamak (MAST) and compared against a previous analysis of helium 
impurity transport on MAST. The impurity density profiles are measured on the low- 
field side of the plasma, therefore this paper focuses on light impurities where the 
impact of poloidal asymmetries on impurity transport is predicted to be negligible. A 
weak screening of carbon and nitrogen is found in the plasma core, whereas the helium 
density profile is peaked over the entire plasma radius. Both carbon and nitrogen 
experience a diffusivity of the order of 10 nr/s and a strong inward convective velocity 
of ~ 40 m/s near the plasma edge, and a region of outward convective velocity at mid¬ 
radius. The measured impurity transport coefficients are consistent with neoclassical 
Banana-Plateau predictions within p < 0.4. Quasi-linear gyrokinetic predictions of the 
carbon and helium particle flux at two flux surfaces, p = 0.6 and p = 0.7, suggest that 
trapped electron modes are responsible for the anomalous impurity transport observed 
in the outer regions of the plasma. The model, combining neoclassical transport 
with quasi-linear turbulence, is shown to provide reasonable estimates of the impurity 
transport coefficients and the impurity charge dependence. 


1. Introduction 

Tokamaks contain a range of partially and fully ionised impurity ions. Helium ash will 
be unavoidable in D-T plasma, while plasma facing materials in ITER will consist of 
tungsten and beryllium. Additionally, seeded elements, such as argon, neon or nitrogen, 
have been identified as suitable mantle and divertor radiators jT] . A build-up of these 
impurities in the main chamber plasma may result in fuel dilution and radiative cooling 
of the plasma core, degrading the tokamak performance. Predicting and controlling 
impurity transport in next-step tokamaks requires the rigorous testing of suitable models 
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against experimental data and an assessment as to how impurity transport depends 
on plasma parameters. Several attempts have already been made in this direction 
(see Guirlet et al. [2] for example), and this paper focuses particularly on the charge 
dependence, in the case of a spherical tokamak (ST). 

The atomic number, Z, and mass number, A & mi/m p , both strongly influence 
neoclassical and turbulent impurity transport. This paper focuses on fully ionised 
impurity ions with Z < 7 and A/Z ~ 2. In this limit, turbulent diffusion can either 
decrease or increase with Z, depending on the character of the dominant turbulent modes 
[3]. The turbulent convective velocity includes contributions from thermodiffusion, 
rotodiffusion and parallel compression [3]. The first two terms decrease and increase 
respectively with Z, while parallel compression is proportional to Z/A and is therefore 
approximately constant for fully ionised impurities. It is worth noting however that 
the magnitude of the turbulent diffusion and convective velocity is expected to saturate 
for fully ionised impurity ions with Z > 10 [3j. In stationary, low collisionality plasma 
neoclassical diffusivity decreases with Z [1] , while in highly rotating plasma neoclassical 
diffusivity becomes a complex function of Z and A [5]. Neoclassical convective velocity 
scales linearly with Z in all collisionality regimes [4j and can therefore dominate over the 
equivalent turbulent rates for highly charged ions, whereas turbulent diffusion typically 
always dominates over neoclassical diffusive rates. 

In light of this complex Z dependence of impurity transport, a large number of 
experimental studies have focused on measuring low-Z (Z < 10) [6H8j, medium -Z 
(10 < Z < 18) P ITT] and high-Z (Z > 18) [T2l ITT] impurity transport coefficients. The 
majority of these experiments have been carried out on conventional tokamaks with only 
a limited number of experimental studies on STs [T51 - H8] . The magnetic field and ion 
temperature is typically lower in STs than in conventional tokamaks, and so ST data 
tests impurity transport models in plasmas where the neoclassical transport is higher and 
more dominant. Additionally, different modes can drive turbulent transport and long 
wavelength ion temperature gradient (ITG) modes, often associated with anomalous 
transport on conventional tokamaks, are frequently stabilised in STs by sheared toroidal 
flows [19] . It was found in a previous study that helium impurity ions are sensitive 
to trapped electron mode (TEM) turbulence in L-mode discharges in the Mega-Amp 
Spherical Tokamak (MAST) jT8] • This paper extends that study using measurements of 
carbon and nitrogen transport made during L-mode MAST plasma scenarios. Nitrogen 
and carbon are analysed in this paper to assess the transport of extrinsic and intrinsic 
impurities respectively, and the Z dependence of impurity transport in these MAST 
plasmas is obtained by comparing with the helium results in fl8] . 

Experimental transport coefficients are determined using the 1.5-D transport code 
SANCO [20] to reproduce the density profile evolution of C 6+ and N' + , measured following 
a short impurity gas puff using charge exchange (CX). Predictions of the carbon, 
nitrogen and helium [18] transport coefficients were obtained using the neoclassical code 
NEO [2DE2] and the gyrokinetic code GKW [23]. The results suggest that, compared 
to helium, carbon and nitrogen have higher rates of inward convective velocity and 
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diffusion near the plasma edge, and the convective velocity of carbon and nitrogen 
(but not of helium) changes sign from inwards to outwards at mid-radius. Carbon, 
nitrogen and helium diffusivities agree within error bars in the core. Our theory based 
impurity transport model suggests that the trend in Z near the edge is caused by TEM 
turbulence, and finds that the reversal in the convective velocity for carbon and nitrogen 
at mid-radius is consistent with neoclassical transport in the Banana-Plateau regime. 

This paper is organised as follows. A description of the L-rnode 900 kA plasma used 
for the analysis is given in section [21 where we also provide details of the diagnostics 
used to measure the CX spectral intensities and the bulk plasma parameters that are 
required as inputs to the simulation codes. In section [3] we present the model of the C 6+ 
and N' + density evolution, and describe the procedure used to obtain the experimental 
transport coefficients. We discuss the measured impurity transport coefficients and 
their dependence on Z in section 01 and interpret these results in terms of an impurity 
transport model based on neoclassical and quasi-linear gyrokinetic theory. Lastly, we 
provide a brief conclusion in section [5] 

2. Plasma Scenario and Diagnostics 

We have chosen to study the dependence on Z of impurity transport in the L-mode 900 
kA MAST discharges described in our initial report [TSJ. This plasma is free from edge 
localised modes (ELM), while TEM turbulence is expected in the outer region of the 
plasma. The transport analysis time window is reduced to a few tens of milliseconds 
during the current flat-top (beginning at f ~ 0.2 s) to avoid the magnetohydrodynamic 
mode (MHD) activity occurring at t ~ 0.28 s. A summary of the key plasma parameters 
are as follows: plasma current I p = 900 kA, toroidal magnetic field Bt = 0.55 T, NBI 
heating power Pnbi = 2.1 MW, on-axis electron density n e = 3.5 x 10 19 m -3 , on-axis 
electron temperature T e = 1.0 keV, I p flat-top and diffusion duration is of the order of 
100 ms, global energy confinement time is of the order of 10 ms, major radius i ? 0 = 0-83 
m, minor radius a = 0.6 m, elongation k = 1.93 and triangularity 6 = 0.4. Plasma 
parameters vary by less than 20 % during the I p flat-top in this particular discharge. 
Time traces of various equilibrium quantities are shown in [18]. 

Helium (He), methane (CH 4 ) and nitrogen (N 2 ) gas were each injected into the 
plasma at the beginning of the I p flat-top using a piezo-valve located on the inboard 


Table 1 . Details of the injected gas species during repeated MAST discharges. 


Mast Shot # 

Injected Gas Species 

Fuelling Duration / ms 

28052 

None 

— 

28053 

ch 4 

35 

29261 

He 

15 

29427 

None 

— 

30439 

n 2 

15 
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Figure 1 . Profiles of (from top to bottom) the electron and ion density, electron 
and ion temperature, ion toroidal velocity and (deuterium) Mach number, effective 
charge, and the safety factor, as a function of the normalised toroidal magnetic flux, 
P = \Z4>n- Each radial point represents the weighted mean over the five repeated shots 
over the time interval 0.21 s < t avg < 0.27 s, with the error bar indicating the standard 
deviation. Both the ion density and effective charge have been calculated using CX 
concentrations of carbon (~ 0.1 %) and helium (~ 1 %). 


side of the lower divertor. The transport analysis time window is taken during the I p 
flat-top time between 0.21 s < t avg < 0.27 s, which is free from the long-lived mode 
and sawtooth MHD activity. The plasma scenario was repeated three times to study 
one impurity gas puff (of the order of 10 ms) per pulse. Two further repeat pulses 
were carried out with no impurity gas injection for comparison. A summary of the 
shot numbers and injected impurity gas is given in table [0 A longer fuelling time was 
needed for CH 4 because the measured fuelling efficiency for carbon was found to be 
lower than that for helium EH The impurity injections can be considered trace since 
they increased the effective charge of the plasma, Z e //, by < 5 %. 

The transport analysis presented in section |4] is supported by simulation codes 
which require, as input, a description of the bulk plasma equilibrium. MAST is equipped 
with a comprehensive set of diagnostics. The Thomson scattering system measures the 
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electron temperature T e and density n e profiles with a radial resolution and sampling 
rate of 10 mm and 240 Hz respectively [25j. Spectrally resolved CX measurements of 
C VI from the CELESTE-3 diagnostic are used to deduce the ion temperature T t and 
toroidal rotation profiles with radial and temporal resolutions of 1 cm and 4 ms [26]. 
Radial profiles of the He 2+ and C 6+ concentration (typically of the order of ~ 1 % and 
~ 0.1 % respectively) are derived from calibrated CX measurements from the RGB 
filtered imaging diagnostic [27] and used to determine self-consistent profiles of Z e ff 
and the main ion density n* 0. Radial and temporal resolutions of 1 mm and 5 ms 
respectively are available from RGB, although typically the radial resolution is reduced 
to 1 — 5 cm in order to improve statistics. The magnetic equilibrium, and hence safety 
factor g, has been reconstructed using the EFITH—b code [28], with D a and motional 
Stark effect (MSE) constraints [29] . 

Radial prohles of several of the plasma parameters described above are shown in 
figure [I] as a function of p = y 7 ^, where (pjy is the normalised toroidal magnetic 
flux. Each radial point represents the weighted mean, with a corresponding standard 
deviation, taken over time ( t avg = 0.21 — 0.27 s with 5 ms resolution), shot (5 repeat 
discharges) and space. Significant variations in the transport coefficients are not evident 
for spatial resolutions of < 5 cm, therefore the plasma parameters are spatially resolved 
to 5 cm. Uncertainties of < 20 % and < 40 % are found for the weighted mean values 
and their radial gradients respectively. The standard deviation of the radial gradient 
was calculated using the method of propagation of errors. 

Measurements of the active CX C VI n = 8 -* 7 (A = 529.1 nm) and N VII 
n = 9 —y 8 (A = 566.9 nm) spectral line intensities, induced by NBI, have been made 
with a spectrometer and multi-chord setup. For carbon we use the existing CELESTE- 
3 diagnostic [26] (primarily measuring T) and V^) consisting of 64 active and passive 
tangential fibres with a spatial resolution of ~ 1 cm. A similar setup is used for nitrogen 
but with 16 active and passive fibres with a spatial resolution of ~ 5 cm. Integration 
times of 5 ms are used for both measurements. 

3. Impurity Density Evolution 

The C VI and N VII CX spectral line intensities are converted into C 6+ and N' + density 
prohles using the CX effective emission coefficients interpolated from the Atomic Data 
Analysis Structure (ADAS) [30] and simulations of the line-integrated neutral beam 
density tat For the impurity transport analysis, only the evolution of the impurity 
density profile shape matters as explained later in this section; absolute concentrations 
are not strictly necessary. Therefore the systematic error from the thermal beam halo, 
which can increase the CX intensity by a constant radial fraction of ~ 20 % for the same 
impurity density [24], will not modify the derived transport coefficients. Secondary 
impurity ‘plume’ emission, induced by electron collisional excitation of the partially 
ionised impurity ‘plume’ ions left over following the CX process, can alter the radial CX 
intensity profile shape. However, we expect any changes in the profile shape caused by 
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Figure 2. Profile evolution of the (a) N' + and (b) C 6+ impurity density profile 
normalised to the maximum value. The start time of the gas puff is subtracted from 
the time axis in both graphs. Both radial flux labels p and r/a are indicated above 
and below the plots respectively. t gas corresponds to (a) 195 ms and (b) 175 ms. 


impurity plume emission to be negligible due to the relatively long ionisation lengths and 
low electron excitation rate coefficients associated with the plume ions [21]. The C 6+ 
and N' + density evolution after the gas puff is illustrated in figure [2j Intrinsic impurity 
density profiles, calculated during the two reference discharges, have been subtracted. 
A region of ~ 20 ms where no signal from the gas puff is measured is found immediately 
after the piezo-valve is triggered. This corresponds to the finite time taken by the 
impurity gas to travel through the plenum and into the main chamber. 

The 1.5-D radial transport code sanco [2D] is used to numerically solve the 
continuity equation for the density of each impurity ionisation stage with a given set 
of boundary conditions, sanco uses a cylindrical geometry, therefore the transport 
equations are written as 


S. 


i,z 


dn IyZ _ Id 
~dT~ -rfr {rTl ’ z) 

Ti,z = ~ + V r np Z 


( 1 ) 

( 2 ) 


where Dj and Vj are the diffusion and convective velocity coefficients for impurity species 
/ and r = \/Vfi ux /2iT 2 R 0 is the effective radius of the flux surface enclosing a volume, 
Vfiux. Ionisation and recombination rate coefficients provided by ADAS are used to 
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Figure 3. A profile of the intrinsic logarithmic C 6+ density gradient from the 
discharges without impurity gas puffs. Each radial point represents the weighted mean 
over time ( t aV g = 0.21 — 0.27 s), space and shot (2 repeat discharges with no impurity 
gas puff), with the error bars denoting a corresponding standard deviation. 


determine the source term, SV,z- Determining the transport coefficients requires a least- 
squares minimisation of the error between the experimental and modelled fully ionised 
impurity density evolution by the appropriate optimisation of Dj and Vi/Dj. 

The SANCO model assumes that the impurity transport coefficients are constant 
in time and equal for each impurity ionisation stage. The first assumption is valid for 
the time window of t avg because the plasma is approximately in steady-state during the 
plasma current flat-top on MAST. In theory, the latter assumption is not strictly valid, 
since our main hypothesis is that there is a charge dependence on impurity transport; 
however this simplification is valid in this case because the fully ionised impurity ions 
are dominant over the majority of the plasma radius. Ionisation timescales for the 
partially ionised impurity ions that exist near the cooler plasma boundary are such that 
any differences in plasma transport between each ionisation stage will have a negligible 
effect on the derived impurity transport coefficients. 

Radial impurity density profiles from SANCO are averaged over flux surfaces and 
hence do not account for any poloidal variation. Low-held side (LFS) localisation of 
impurities is an important issue for transport studies on STs because a low aspect ratio 
and strong toroidal rotation, typical of ST plasma, can enhance the poloidal asymmetry 
that arises from the outward centrifugal force [31]. Unfortunately, CX impurity line 
intensities fall below the detection threshold of the diagnostics on the high-held side 
(HFS) of the plasma due to the exponential fall-off in neutral beam density, therefore 
the poloidal variation of the impurity density cannot be measured. However, a prediction 

















Figure 4. A comparison between the temporal evolution of the experimental and 
simulated SANCO C 6+ density profile for a number of different radii. The blue shaded 
region indicates the timing of the impurity gas puff. 


of the LFS localisation can be made using [3Tj 

^-r = (l - ^ AA) (3, 

< n I>z > rrii \ m z T e + J 

where n FF z s is the impurity density at the low-held side of the plasma, e is the inverse 
aspect ratio and < n I)Z > is the average impurity density on a flux surface. With 
profiles of Mi, Z eff, T e and T % from figure [1] values of Uq FS / < nc6 ><0.1 are found 
and therefore any change in the derived carbon transport coefficients is expected to be 
within the fitting error. The poloidal variation of heavier impurities in similar plasma 
conditions may be significant; for example, a trace injection of argon with dominant 
ionisation state Ar 16+ could have values of n\ F f 6 / < riArw >~ 0.4. 

To fit the transport coefficients, first we parameterise Di and Vj/Dj over 10 equally 
spaced radial points corresponding to a spatial resolution of 5 cm. Fixed values of 
Vi/Di are used for the first iteration of the minimisation with Dj as a free parameter. 
Measurements of the C 6+ logarithmic radial density gradient made during the two 
reference discharges provide an estimate of the zero-flux (T j >z = 0) peaking factor 
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and hence an initial guess for Vj/Dj: 

1 d In rij.z 


Vj/Dj = 


n I,Z 


dr 


(4) 


The measured zero-flux peaking factor profile shown in figure [3] suggests an inward 
pinch in the region 0.5 < p < 1.0 and a region of outward convective velocity between 
0.3 < p < 0.5. Next we set both Dj and Vj/Dj as free parameters allowing for any 
correction to Vj/Dj. After Dj and Vj have been determined over the equally spaced 
radial grid, we then remove grid points within the regions of approximately constant D] 
and Vj. The errors, illustrated by the shaded region in figure El represent one standard 
deviation for Dj at constant Vi/Dj and for Vi/Dj at constant Dj. 

The edge model in SANCO requires inputs for the neutral influx T 9as , the recycling 
coefficient R rcy , the parallel loss time 7 ]|, the neutral kinetic energy E th , the fuelling 
efficiency €f and the values of Dj OL and Vf OL in the scrape-off-layer (SOL) pTJ. T 9as 
is calculated by integrating the measured impurity density over the plasma volume. 
We set R rcy = 0.01 to represent the low recycling typical of carbon and t\\ = 1 ms 
which was based on the SOL connection length and Eth- The value of E t h determines 
the distance that the neutral particles travel into the plasma before ionising: e.g. with 
E t h = 2 eV the neutrals travel ~ 10 cm into the plasma (corresponding to region between 
0.8 < p < 1.0). Underestimating E th will cause Dj and Vj to be overestimated over the 
range determined by the maximum penetration depth of the neutrals. 

Atoms and molecules leaving the gas valve have thermal energies corresponding 
to room temperature. As they interact with the scrape-off-layer (SOL) plasma, the 
molecules dissociate at a threshold T e forming atoms with Franck-Condon breakup 
energies of a few eV. The atoms are kinetically energised by elastic collisions and CX 
reactions with plasma ions, however their ionisation times are too short for them to 
thermalize fully with the SOL plasma (generally around 10 - 20 eV on MAST) and 
therefore we assume that E th ~ 2 eV [H1I32]. We note that increasing E th to around 
10 eV caused a reduction of Vj and Di within the region 0.7 < p < 1.0, but this value 
of E t h is much higher than that in MAST. The values of ep, Dj OL and Vf OL were 
tailored to match the magnitude of the SANCO and CX impurity density profiles up to 
the expected penetration depth of the neutral impurity atoms ( 0.8 < p < 1 . 0 ). 

A comparison of the SANCO and experimental impurity density evolution is shown 
in figure [H The impurity particle flux perturbation induced by the gas puff was not 
evident within the plasma core p < 0.2, therefore the transport analysis is restricted to 
the spatial range of 0.2 < p < 0.8. As a goodness of fit metric, we use the reduced chi- 
squared statistic, Xr ed = (/* — Vi) 2 1 o 2 j DOF where i is the number of measurements, 

a is the variance and DOF is the degrees of freedom. The C 6+ density fit (between 
0.21 < t < 0.27 s) produced a x 2 e d = 5.2, while the N 7+ density fit was quantified 
by Xred = 4-4. A model in accordance with the variance should give x 2 e d = 1) while 
Xred >> 1 indicates that the model is incorrect and Xred > 1 indicates that a has been 
underestimated. For both C 6+ and N 7+ we set a = 0.2 however, assuming DOF cs 100, 
a Xred = 5-2 would suggest a more realistic (but still modest) variance of ~ 2 cr. 
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Figure 5. Radial profiles of impurity transport coefficients for helium, carbon 
and nitrogen: (a) the diffusivity, D- (b) the convective velocity, V; and (c) the 
corresponding steady state impurity peaking factor, —V/D. NEO simulations of the 
neoclassical impurity transport coefficients for helium and carbon are shown as dashed 
and dashed-dotted lines. The quasi-linear and neoclassical transport coefficients for 
helium and carbon have been added together at p = 0.6 and p = 0.7 as shown by the 
symbols. The symbols within p < 0.5 indicate the typical neoclassical error. 


4. Charge Dependence 

The Dj, Vj and Vi/Dj sanco profiles are shown in figure [51 Carbon and nitrogen 
transport coefficients agree within error bars over the analysed radial region. Both 
experience a strong inward pinch (0.5 < p < 0.8), a region of weak outward convective 
velocity around mid-radius (0.4 < p < 0.5), and a diffusive region in the plasma core 
(p < 0.4). Compared to helium, they experience a larger inward pinch and diffusivity 
near the plasma edge and a change of convective velocity direction (from inwards to 
outwards) at mid-radius. This suggests that higher Z impurities in an ST will stagnate 
around the plasma mid-radius, whereas lower Z impurities will accumulate further into 
the plasma core. 
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We model the neoclassical particle flux at each radial point using the local drift 
kinetic neoclassical code neo [21][22]. Toroidal rotation of the bulk plasma and self- 
consistent concentrations of C 6+ and He 2+ from CX measurements are included in the 
NEO inputs. Neoclassical predictions for helium and carbon are shown by the dashed 
and dashed-dot lines respectively in figure El rates for nitrogen differ little from carbon 
and are therefore omitted. In the region 0.4 < p < 0.7, neoclassical convective velocity 
and diffusion are lower in magnitude than experiment for both impurities. Carbon 
and helium diffusivities agree with the neoclassical predictions in the region p < 0.4. 
Neoclassical simulations reproduce the reversal of the sign of the impurity convective 
velocity with radius, which is observed for carbon but not for helium. At p ~ 0.4 both 
C 6+ and He 2+ ions are in the low collisionality Banana-Plateau (BP) regime and hence, 
assuming that X) = T/, the neoclassical BP convective velocity is (see Eq. 3 of [4j) 


ybp 


d) p z 


d In n t 
dr 


+ 



1 

~Z 


1 


1 ^ 1,2 M ,2 j & In Ti 

Z k[ , k l xl J dr 


(5) 


where D bp is the BP diffusivity and kf ,j are the viscosity coefficients. The direction of the 
convective velocity is governed by the density and temperature logarithmic gradients, 
and by the factor in brackets multiplying the temperature gradient. To illustrate the 
dependence on Z of this multiplying factor, we use the asymptotic form of kfj (defined 
in Eq. A. 15 of [3]) with the dominant contribution from the banana regime. In this 
limit, the multiplying term reduces to ~ —0.19 and ~ —0.34 for helium and carbon 
respectively. Since the temperature gradient is also negative, increasing Z therefore 
enhances the screening of the impurity ions providing an explanation for the reversal in 
the impurity convective velocity direction observed for carbon. 

The turbulent particle flux is modelled at p = [0.6, 0.7] with the local flux-tube 
gyrokinetic code gkw [23] 0. These two flux surfaces are in the steep density gradient 
region where TEMs are linearly unstable [18]. Turbulent contributions within p < 0.5 
have not been modelled since the flatter density gradient stabilises the TEMs. A pure 
plasma (n e = rii) has been assumed for the GKW simulations since the low impurity 
concentrations ( Z e ff < 1.3) have negligible effects on the quasi-linear particle fluxes [ 33] , 
Linear growth rate 7 spectra calculated with gkw for p = [0.6, 0.7] are shown in the 
top panels of figure El as a function of k y pi, where k y is the binormal perpendicular 
wavenumber and pi is the local Larmor radius of the main ion. The stabilising effect 
of equilibrium toroidal flow shear = u'B e /B ~ u'e/q , where u' is the radial gradient 
of the equilibrium toroidal rotation has been determined with the gyrokinetic code 
GS2 [34|[35] using the technique described in [THUSBjI - The bottom two panels of figure 
[6] show the mixing length diffusivity estimate 7/A; 2 and the mode frequency c 0 . The 
dominant modes, corresponding to the wavenumber of the maximum value of 7/A; 2 , are 
within the TEM region; k y pi ~ 3.0 and k y pi ~ 2.0 for p = 0.6 and p = 0.7 respectively. 

t Input files used for the NEO and gkw simulations can be found at 
https://bitbucket.org/gkw/publication-input/src/HEAD/2015_Henderson_PPCF/ 
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Figure 6. Linear growth rate spectra (top panels), mixing length diffusivity estimates 
(middle panels) and mode frequencies (bottom panels) are shown for (a) p = 0.6 and 
(b) p = 0.7. Calculations were performed with and without flow shear using the GS2 
(triangles) and GKW (diamonds) codes respectively. The dashed lines in the top panels 
represent the value of the flow shear (normalised to c s /a), while the shaded region 
in the middle panels represents the effective electron heat diffusivity calculated using 
JETTO [371138] . 


We choose the value of k y pi = 2.6 for both flux-surfaces to use in our quasi-linear analysis 
described below. 

The modelled particle flux can be written as a linear summation of the neoclassical 
and turbulent contributions, 


•HMODEL 

1 I,Z 


7WNEO T^turb 

1 -|- L _b Vr NEO + Vr UTb 


ni,Z 


( 6 ) 


where 1/L a = — d In a/dr. The diffusive and convective terms are extracted from the 
simulation output hies by including trace species with different gradients in the input 
files. In the quasi-linear approach we model the saturated turbulence as the dominant 
linear mode with an amplitude that reproduces the effective electron heat diffusivity, 
X e J ^, calculated using the interpretative transport code JETTO [37]l38] (see table [2]). 
This allows us to determine the turbulent impurity particle diffusivity and pinch from 


their corresponding GKW values, D GKW and E GKW , and electron heat diffusivity, x GKW , 111 
linear GKW simulations as: 


j-^turb _ 


yGKW 

Ae 


D c 


yturb _ 


yGKW 

Ae 


v c 


(7) 


At steady state (i.e. 
peaking factor: 


with T i z = 0) Eq. O provides the modelled dimensionless 


turb 


J ni,z 


= —a— - —t (8) 

£)NEO _|_ jjturb > 

where we choose the value of D GKW , V GKW and Xe KW associated with the dominant 
linear mode. Strictly, the quasi-linear single value of each of these profiles must be 
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Figure 7. Quasi-linear carbon and helium transport coefficients from GKW as a 
function of normalised wavenumber for two flux surfaces, p = 0.6 and p = 0.7. The 
shaded region represents the dominant linear mode numbers for each surface. 


associated with the same k y pi. In most conventional tokamak discharges, the ITG modes 
responsible for impurity particle transport also drive the anomalous heat diffusivity and 
hence the above method is valid and used in many studies [39l ITT] . However on MAST 
it is the electron turbulence that is thought to drive the anomalous heat transport in 
most discharges [36]. We therefore assess the applicability of this quasi-linear approach 
based on its agreement with experiment; a rigorous test would include a comparison 
with non-linear gyrokinetic simulations which are not attempted here. 

Figure [7] shows the wavenumber profiles of quasi-linear and V} urb for radial 

locations p = 0.6 and p = 0.7. The impurity transport coefficients peak around k y pj ~ 1 
(noting that k y pc ~ k y pi/V2 and k y pue ~ k y pi/\/6 ). Turbulent diffusion is moderately 
higher for carbon compared to helium; a trend previously found using a reduced fluid 
model for TEMs [3]. In this fluid description, the diffusivity depends upon Z directly 
and indirectly due to the fluid curvature. For impurities in the trace limit, the turbulent 
convective velocity can be decomposed as [42] 


-yturb 


x 


eff nGKW 

U I 


yGKW 

Ae 


C T — + C u v! + C p 

L/t 


(9) 


where Ct , C u and C p are the thermodiffusion, rotodiffusion and convective dimensionless 
turbulent transport coefficients respectively [3]. V turb moderately increases with Z 
because C u is dependent on the impurity Mach number; the ratio of the two impurity 
Mach numbers is approximately Mc/M He « yj me / mHe ~ 1.7. We also noted that 
Ct was moderately smaller for carbon compared to helium, however this decrease was 
modest compared to the increase in C u . There is little change in C p since this depends 
upon the ratio Z/A, which is constant for fully ionised impurities. 

To determine the uncertainties in the modelled transport coefficients, the 
simulations have been repeated whilst systematically varying each experimental input 
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according to their standard deviation. The biggest uncertainty in Zl/ E0 was caused by 
the q profile (A q/q ~ 15 %), which is unsurprising considering that both low and high 
collisionality neoclassical diffusivities are proportional to q 2 . Vf E0 was equally sensitive 
to the error in the temperature and density radial gradient (~ 40 %). Bremsstrahlung 
measurements give Z e ff profiles that are typically higher than from CX measurements 
in similar L-mode plasmas, with Z e ff increasing up to ~ 2 in the region p > 0.5 |2T] - We 
explored this uncertainty by increasing the CX C 6+ concentration by a factor of 10 in the 
plasma edge and redoing the NEO simulations, but the corresponding errors were modest 
compared to those arising from uncertainties in q and in the density and temperature 
gradients. The main source of uncertainty in the gyrokinetic diffusivity arises from the 
error in x e e (~ 30 %), while the convective velocity uncertainty is sensitive to both x^J 
and the ion temperature gradient (which multiplies the thermo diffusion coefficient). 

Modelled transport coefficients for helium and carbon at p = [0.6, 0.7], including 


Table 2. Comparison of the experimental and modelled transport coefficients for 
carbon and helium at p = [0.6, 0.7]. The individual contributions to the model from 
NEO and turbulence are given in the last two columns of the table. 



Experiment 

Model 

Total 

NEO 

Turb 

p = 0.6 

Xe / m 2 /s 

4.4 ± 1.3 

— 

— 

— 

D c / m 2 /s 

1.8 ± 0.3 

2.1 ± 0.5 

0.6 

1.5 

D He / m 2 /s 

1.5 ± 0.5 

1.4 ± 0.2 

0.7 

0.7 

D c /D He 

1.2 ± 0.4 

1.5 

0.8 

2.1 

V c / m/s 

-12.5 ± 2.4 

-10.5 ± 5.7 

0.1 

-10.6 

V He / m / s 

-5.5 ± 1.2 

-7.4 ± 3.3 

-1.0 

-6.4 

V c /V He 

2.3 ± 0.7 

1.4 

-0.1 

1.6 

—aVc/D c 

4.0 ± 0.5 

3.0 ± 1.8 

-0.1 

4.2 

—aVne/D He 

2.2 ± 0.7 

3.2 ± 1.5 

0.9 

5.3 

p = 0.7 

Xe / m 2 /s 

6.1 ± 2.1 

— 

— 

— 

D c / m 2 /s 

4.9 ± 0.9 

3.1 ± 0.7 

1.0 

2.0 

D He / m 2 /s 

2.2 ± 1.0 

2.5 ± 0.4 

1.4 

1.1 

Dc/DHe 

2.3 ± 1.1 

1.2 

0.8 

1.8 

V c / m/s 

-34.2 ± 9.7 

-13.5 ± 5.4 

-3.9 

-9.6 

V He / m/s 

-12.7 ± 4.9 

-9.3 ± 3.6 

-2.3 

-7.1 

V c /V He 

2.7 ± 1.3 

1.4 

1.7 

1.4 

—aVc/D c 

4.1 ± 0.3 

2.7 ± 1.2 

2.3 

2.9 

-aV He /D He 

3.3 ± 1.4 

2.3 ± 0.9 

1.0 

3.8 






15 


both neoclassical and turbulent contributions, are shown by the square and circle 
symbols respectively in figure [5] and summarised in table [2j Note that symbols are also 
shown (for carbon) within p < 0.5 to indicate the typical error bar of the neoclassical 
impurity transport coefficients. The magnitude of the carbon and helium transport 
coefficients has been reproduced within the uncertainties for p = 0.6, however the 
magnitude of the experimental convective velocity for carbon is significantly higher 
than the model prediction at p — 0.7. Encouragingly, the modelled trend in Z for 
both diffusion and convective velocity rates, given by + -D£ KW )/(_Dj/ E ° + Df]™) 

and similarly for the convective velocity, agrees within the experimental uncertainty of 
D c /D He and Vc/Vue at both p = 0.6 and p = 0.7. Note that the uncertainty in the 
ratio of the modelled diffusivity for carbon and helium (or in the corresponding ratio 
for the modelled convective velocities) is small because of error cancellation. 

Most noticeable from the experimental impurity transport coefficients shown in 
figure [5l is that the dependence on Z agrees with the trends expected for TEM 
turbulence that is predicted to be unstable. This is opposite to the trend expected 
for ITG turbulence, where turbulent diffusion is expected to decrease with Z, while 
both thermodiffusion and rotodiffusion are directed outwards [31143]. Thus, the main 
conclusion from this work is that a combination of neoclassical and TEM driven 
transport is most likely causing the observed impurity particle flux in the outer region 
of the plasma. Furthermore the model, combining neoclassical transport with quasi- 
linear turbulence, is shown to provide reasonable estimates of the impurity transport 
coefficients over the analysed spatial region. 

5. Conclusion 

This paper has evaluated fully ionised carbon and nitrogen impurity transport 
coefficients during two repeat L-mode 900 kA MAST discharges for comparison with 
earlier results for helium. There is little difference between the carbon (intrinsic) and 
nitrogen (extrinsic) transport coefficients which is not surprising, considering that both 
species are fully ionised over the majority of the plasma radius and therefore have a 
similar atomic charge and mass. Both experience a diffusivity of the order of 1-10 m 2 /s 
and a strong inward convective velocity of ~ 40 m/s near the plasma edge, and a region 
of outward convective velocity at mid-radius. In comparison with helium, carbon has 
higher diffusivity and convective velocity near the plasma edge and furthermore its 
convective velocity has a reversal at mid-radius. 

Neoclassical and quasi-linear gyrokinetic simulations of the particle flux have been 
run with NEO and GKW respectively. Neoclassical transport alone is sufficient to explain 
the observed impurity transport of each species within p < 0.4, but cannot explain the 
magnitudes of the transport coefficients or trend in Z in the region 0.4 < p < 0.8. 
The reversal of the convective velocity direction at mid-radius is consistent with 
neoclassical Banana-Plateau transport. A quasi-linear approach has been used to 
calculate the magnitude of turbulent diffusivity and convective velocity, and indicates 
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that in these MAST discharges the dominant linear mode is in the trapped electron mode 
wavenumber region (because of the suppression of ITG turbulence by sheared toroidal 
flow). Summing the neoclassical and quasi-linear turbulent contributions to the particle 
flux provides a model of the impurity transport that reasonably describes its magnitude 
and dependence on charge; however it is noted that the observed carbon transport 
at p = 0.7 is still significantly higher than the model. For helium, the anomalous 
convective velocity is dominated by thermodiffusion, whereas both rotodiffusion and 
thermodiffusion contribute to the anomalous carbon convective velocity. In summary, 
this paper provides further support for the existence of trapped electron modes in MAST 
plasmas and their impact on low-Z impurity transport. 
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